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SUMMARY 

An approach for developing exact one-dimensional 
conduction-convection finite elements is presented. Exact 
Interpolation functions are derived based on solutions to the 
governing differential equations by employing a nodeless 
parameter. Exact interpolation functions are presented for 
combined heat transfer In several solids of different shapes, 
and for combined heat transfer In a flow passage. Numerical 
results demonstrate the exact one dimensional elements offer 
significant advantages over elements based on approximate 
interpolation functions. 

1, INTRODUCTION 

A significant development in finite element methodology 
for forced convection analysis is the concept of upwind 
weighting functions [1]. For the one-dimensional convective-, 
diffusion equation upwind weighting functions are expressed 
In terms of an upwind parameter which can be varied to adjust 
the accuracy of the finite element solution. Optimal values 
of the upwind parameter lead to exact values of the tempera- 
ture at nodal points. Other elements in one- dimensional 
thermal and structural analysis also predict exact nodal 
values. For example, one dimensional solutions [2] for heat 
conduction with internal heat generation by linear interpola- 
tion functions yield exact temperatures at the nodes. The 
ability of finite elements to predict exact nodal temperatures 
has previously -been regarded as a property of a few particular 
equations and to be of limited generality. There is, however, 
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a broad class of' one'dimensional problems governed by linear 
ordinary differential eguations for which finite elements can 
be developed to yield exact solutions. The one-dimensional 
convective-diffusion equation and the one-dimensional conduc- 
tion equation are special cases of this class of problems. In 
a previous paper [3]. the authors have presented two examples 
of one-dimensional thermal -structural finite elements which 
yield exact solutions for linear problems. 

The purpose of the paper is to present an approach for 
developing one-dimensional conduction-convection finite ele- 
ments which yield exact nodal temperatures and an exact 
variatior of temperature within an elenenc for steady-state 
linear analysis. The one-dimensional conduction-convection 
problems considered are first discussed. A nodeless variable 
approach for deriving exact interpolation functions is then 
presented and illustrated for several cases. Next, the gener- 
al form of the element equations is presented and discussed. 
Finally, the benefits of the exact finite element formulation 
are illustrated for two problems by comparing numerical solu- 
tions from finite elements with exact and linear interpolation 
function?. 

2. EXACT ELEMENT FORMULATION 

2.1 Governing Equations and Boundary Conditions 

The geometry and terminology for eight one-dimensional 
conduction-convection cases are shown in Fig. 1. Cases 1-7 
are solids of various shapes where heat transfer consists of 
conduction which may be combined with: (a) surface convection, 
(b) internal heat generation, and (c) surface heating. Case 
(8) is a one-dimensional flow where the heat. transfer consists 
of fluid conduction and mass- transport convection which may be 
combined with: (a) surface convection, and (c) surface heat- 

ing. In the figure, 0 is the volumetric heat generation rate, 
q is the surface heating rate, h is the convective heat 
transfer coefficient and T^ is the environmental temperature 
for the convection heat exchange. Heat transfer in Cases 1-8 
is described by differei tial equations of the form 

where a^ , 1 » 0,1, 2, 3, are functions which depend on the 
geometry and thermal parameters of each case. Finite element 
equations for eq. (1) may be formulated by the method of 
weighted residuals, but due to the presence of the odd-order 
derivative an unsynmetrical coefficient matrix will result. 
However, if eq. U) is first cast into self-adjoint form, the 
differential operator will be symmetric. Then, finite element 
matrices can be derived by the method of weighted residuals or 
a variational method, and coefficient matrices will be 


tl. 

Case t ROD Case 2. SLAB 



Case 3. HOLLOW CYLINDER Case 4. HOLLOW SPHERE 



Cose Z SPHERICAL SHELL Case 8. FLOW PASSAGE 


Figure 1 

ONE DIMENSIONAL CONDUCTION AND CONVECTION CASES 


symmetrical, [2]. Eq. (1) Is written in self-adjoint form by 
multiplying it by the factor 


such that 


P • exp[;(a^/aQ) dx] 

• R(x) 


( 2 ) 


where Q(x) * ^2 * *3 


The governing self-adjoint differential equations for 
Cases 1-8 are shown in Table 1. Additional terms not previous- 
ly defined are k the thermal conductivity, p the perimeter of 
the rod or flow passage, A the cross-sectional area of the rod 
or flow passage, m the mas^ flow rate and c the specific heat. 
Table 1 shows the oceffinercs of eg. (2) for each case. Note 
that for each case a diff«tent differential equation may 
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where P * exp(-tncx/kA) 


result by combining conduction with each of the three heat 
loads or combinations of these heat loads. 

The boundary conditions for e<. (2) consist of specifying 
the temperature or the temperature gradient at the endpoints 
of the solution domain, a<x<b. i'Or cases I>8, the boundary 
conditions considered herein are: 

T ■ constant, or 

-k ■ q , or (3) 

-k f • h(T-T^) 


2.2 I nterpolation Functions 

Exact interpolation functions for finite element formula- 
tions of eq. (2) can be derived from the general solution to 
the differential equation. The general solution has the form, 

T(x) • Cj^ fj(x) C 2 f 2 (x) + 6(x) (4) 

where f.(x) and f-(x) are linearly independent solutions of 
the homogeneous equation.c. and C 2 are constants of integra- 
tion, and G(x) is a particOlar soTution. A finite element 
with two nodes is formulated by imposing the conditions 

T(xj^) • Tj T(x2) * T 2 (5) 

where x^ , i ■ 1,2 are nodal coordinates, and T., i • 1,2 are 
the nodil temperature values, To accommodate the particular 
solution Gu)> the element temperature interpolation is 
written as 

T(x) . Nq(x) T^ + N^(x) T^ + N 2 (x) T 2 (6) 


where T is a nodeless parameter. The nodeless parameter T 
is similar to a nodeless variable [2] except the nodeless 
parameter is uniquely determined for each differential equa- 
tion. Note that N (x) is zero at the nodes to satisfy (5). 

The interpolation runctions N . , i -0,1,2 are determined by 
imposing the boundary conditions, eq. (5), on eq. (4) and then 
writing the result in the form of eq. (6). The nodeless para- 
meter is identified by writing the particular solution 
G(x) “ T^ g(x). Thus, in general 

° f2{x.)g(x2)-f2(x2)g(x.) 

NJx) - g(x) * 


^l(x2)g(xj)-fl(x^)g(x2)^ 

[ g ] f. 


(x) 


(7a) 


(7b) 


N^(x) 

N2(x) 

where 


^2(x2)^l(x) * f^(x2)^2^^^ 
B 

f^(x^)f2(x) - f2(x^)fj(x) 

0 * f^(x^)f2(x2) - f^(x2)f2(^l)* 


(7c) 


Element interpolation functions derived in the form of eq. (2) 
are not generally the same as used in conventional finite 
elements since they result from the differential equation 
solution. Interpolation functions, in fact, may range from 
simple polynomials to higher transcendental functions depend- 
ing on the equation. The exact interpolation functions, eqs. 
(7), are derived for specified temperature boundary conditions, 
but they also yield exact solutions for the gradient boundary 
conditions, eq. (3), when these boundary conditions are con- 
sistently incorporated in the finite element equations. Since 
the interpolation functions are exact, temperature gradients 
and fluxes computed from these functions are also exact. 

Nodeless parameters and element interpolation functions 
have been derived for Cases 1-8 with heat loads (a) - (c) , but 
for brevity only typical results are presented herein. Typical 
nodeless parameters appear in TableZand typical interpolation 
functions appear in Table 3. 
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where w « ln(b/a) 


TABLE 3 

ELEMENT INTERPOUTION FUNCTIONS 
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where 


ra * (hp/kA)* 
w * 1n(b/a) 
a » mc/2kA 

6 » (iii^ + a^)^ 



2.3 Element Matrices 


For general interpolation functions, element matrices for 
eq. (2) are 



(8a) 



where K. . Is a typical term in the symmetric conductance 
matrix,^ and F. Is ?. typical term in the heat load vector. For 
a node on the^ boundary additional terms are required to repre- 
sent surface heating or a convective heat exchange. For the 
exact Interpolation functions, element matrices have the 
-general form 
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since 1 » 0,1,2. Because T is known the first equation may 
be uncoupled from the nodal ^unknowns in the second and third 
equations. Thus, the exact element matrices have the same 
size as a conventional linear element and can be written as 
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(9) 


Eq. (9) can be simplified further by noting that K.q-K-q-O 
self-adjoint equations. The last result can be sh6wn oy 
observing that N . , 1*1,2 are solutions of the ordinary 
differential equitlon, (2). Multiplying this equation by N 
and Integrating by parts shows that K.q-ICq* 0. Thus, for 
self-adjoint equations, element matrices need to be evaluated 
only for 1*1,2. 


As emphasized previously, the development of the element 
matrices, eqs. (8), and element equations, eq. (9), are based 
upon the differential equation being written In the self- 
adjoint form, eq. (2). As an alternate approach, exact finite 
element matrices may also be derived from the standard form of 
the differential equation, eq. (1). However, element matrices 



win not be synmetrlcal due to the odd-order derivative, and 
the exact element equations will have K,, 1-1.2. 

The principal advantage of the exact one-dimensional 
finite elements is the superior accuracy In comparison to 
elements based on approximate interpolation functions. The 
principal disadvantage of the exact elements Is the additional 
effort required to form the more complex Interpolation 
functions and to evaluate the element matrices. This disad- 
vantage has been overcome. In part, by using a computer-based 
symbolic manipulation language MACSYHA to perform the algebra 
and calculus required for these derivations. 

3. NUMERICAL EXAMPLES 

Two numerical examples are presented to Illustrate the 
benefits of the exact elements; other examples for Cases 1 and 
3 appear In [4]. The solution of problems with exact elements 
typically requires only a few elements because mesh refinement 
Is not required In regions of large temperature gradients. 

The number of elements used Is determined by changes In 
geometry, material properties or heating variations. After 
nodal temperatures are computed, temperatures are computed at 
several points within each element by using the element Inter- 
polation functions, eq. (6). 

3.1 Coffee Spoon with Conduction and Convection 

Rod elements with conduction and convection (Table 1, 

Case 1(a)) are used to model one-dlmenslon^l heat transfer In 
a coffee spoon (Fig, 2). The lower one-half of the spoon Is 
convectively heated by the coffee at a specified temperature 
of 150*F, and the upper one-half of the spoon Is convectively 
cooled by the atmosphere at a specified temperature of 50®F. 
The ends of the one-dimensional spoon model are assumed to 
have negligible heat transfer. Element matrices for the evxct 
finite elements appear In reference [4]. Temperatures for- the 
spoon model are computed from; (1) two exact finite elements, 
(2) two elements with linear temperature Interpolation, (3) ten 
elements with linear temperature Interpolation. 

The temperature distributions (Fig. 2) show the exact 
temperature distribution computed by the two exact elements 
with three nodal unknown temperatures compared with predictions 
from the elements with linear temperature Interpolation. Two 
linear elements predict the temperatures at the three nodes 
with fair accuracy, but the linear Interpolation functions are 
not capable of representing the zero temperature gradient 
boundary conditions. Temperatures computed from ten linear 
elements, however, show excel! tnt agreement with the exact 
finite element solution. Solutions for similar problems havo 
also shown this trend; typically about five to ten elmnents 
with linear Interpolation per exact element are required for 
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Figure 2 

ROD ELEMENT ANALYSIS OF COFFEE SPOON 
WITH CONDUCTION AND CONVECTION 

3.2 Heat Transfer In Merging Flows 

In merging flows, significant temperature gradients may 
occur at the flow confluence; therefore, merging flows provide 
a good test for finite element evaluations. In reference [4], 
the authors compared conventional (Bubnov^Galerkin) finite 
element solutions with upwind (Petrov-Galerkin) solutions for 
steady-state and transient merging flows. Fig. 3 presents 
the geometry and terminology of a merging flow where conduction 
is combined with mass-transport and surface convection. 

Temperatures for the merging flow (Fig. 3) were computed 
from: (1) three exact finite elements with four nodes, ar:i 

(2) 75 conventional elements with 76 nodes. There is good 
agreement between the two solutions except in the vicinity 
of the flow confluence where the conventional element shows 
oscillations indicating a need for mesh refinement in the 
region of large temperature gradients. The clear superiority 
of the exact elements is demonstrated. Exact nodal tempera- 
tures are also predicted in problems of this type (without 
surface convection) by optimum upwind elements, but a larger 


number of nodes ere re<;u1red since upwind elements do not 
predict exact temperatures within an element. 



Figure 3 

EXACT AND CONVENTIONAL RNITE ELEMENT ANALYSIS 
OF MERGING FLOW WITH CONDUCTION. MASS TRANSPORT 
CONVECTION AND SURFACE CONVECTION. 


4. CONCLUDING REMARKS 

An approach for developing one-*d1ir.ens1onal , exact finite 
elements for linear, steady conduction-convection analysis Is 
presented. The element Interpolation functions employ the 
solution from the governing second-order differential equation by 
utilizing a nodeless parameter approach. The elements require 
two nodes and yield exact nodal temperatures and an exact 
variation of the temperaturewl thin an element. Exact nodal 
Interpolation functions are presented for several cases of 
conduction ard convection. Numerical results are presented 
for conduction and surface convection In a rod element model, 
and for conduction, mass-transport convention and surface 
convection in a merging, one-dimensional flow model. The 
examples demonstrate that for one-dimensional analysis, the 
exact elanents offer significant advantages over conventional 
elements based on approximate Interpolation functions. 
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